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Quasiclassical Green's functions of a short superconductor/normal metal/superconductor 
(SNS) junction are calculated for multi-mode junction with arbitrary electron transmission 
properties. They provide the basis for description of all of the junction transport characteris- 
tics, and are used to calculate the equilibrium spectral density of supercurrent fluctuations. 



1 Introduction 

The process of multiple Andreev reflections (MAR) is the dominant mechanism of electron 
transport in Josephson junctions with large electron transparency. One of the main qualitative 
features of MAR is the coherent transfer of charge by large quanta, the size of which increases 
with decreasing bias voltage V. At low bias voltages, eV <C A, where A is the superconducting 
energy gap, this leads to strong increase of non-equilibrium current noise in the MAR regime, the 
fact that oresently attracts interest to the noise properties of high-transparency junction both 
in theoryllli'i and experiment Si'i. At V smaller than the quasiparticle energy relaxation rate 7, 
eV <C hj, the MAR-noise saturates and turns into equilibrium supercurrent noise that is best 
described in terms of the thermal fluctuations of the occupation factors of the Andreev states 
carrying supercurrent. Since the dynamics of these occupation factors is driven by inelastic 
transitions, comprehensive theory of the current noise in high-transparency junctions should be 
based on the Green's functions formalism which provides a rigorous description of the inelastic 
processes. Until now, such a theory was developed only for the single-mode ballistic junctions^. 
The aim of this work is to calculate quasiclassical Green's functions of a short multi-mode 
Josephson junction with arbitrary electron transmission properties. They are used then to derive 
the general expression for equilibrium supercurrent noise and to find this noise in disordered SNS 
junctions with diffusive electron transport. 

2 Junction model and Green's functions 

The basic model of a short Josephson junction with large electron transparency is a constriction 
between two superconductors that supports N propagating electron modes. If the constriction 
length d is much smaller than the superconductor coherence length ^, the only characteristic of 
the constriction relevant for electron transport is the elastic scattering matrix S: 

;,). (1) 

which is energy-independent in the relevant energy range set by the superconducting energy gap. 
Unitarity of S implies that the N x N transmission and reflection matrices can be expressed as: 
t = u\\'1Du2 , t' = u^y/Du^ , r = u[\'liu4 , r' = —u\^^/Ru2 , where D is the diagonal matrix of 
N transmission probabilities, R = 1 — D, and Uj are unitary matrices. 

Transport characteristics of such junctions can be expressed in terms of their quasiclassical 
Green's functions. Calculation of these functions for the multi-mode junctions presented here 
follows closely similar calculation^ in the single-mode case. The quasiclassical Green's functions 
can be written as 

G{z, z') = Y,{E±e^'^^^'-''^ + F±e±^'=^(^+^')) , (2) 
± 



where z, z' are the coordinates along the current flow, and kp is the Fermi momentum. The 
amphtudes E and F are slowly varying functions of coordinates in the junction electrodes and 
have matrix structure in the transverse mode space. Because of the electron scattering, they 
change rapidly in the junction region ~ 0). Since the junction is short, d <ti ^, rapid variation 
of the amplitudes E and F across it can be described as the boundary condition relating them 
at z ^ —0 and z +0. This condition is derived from the fact that in the junction region, 
G{z, z') satisfy regular Schrodinger equation as a function of z and z', and therefore can be 
expressed in terms of the products of the two scattering states, 

G{z,z')= J2 ^.(^)A.,V,(/). (3) 
Here ipi is the solution of the Schrodinger equation for particles incident from the left, 

and V'2 is similar solution for scattering from the right. Note that "infinity" in eq. @) means 
I z |S> d and still corresponds to z ~ 0, i.e. | z |<C ^, on the length scale of variations of the 
amplitudes E and F in the electrodes that is given by the coherence length ^. Comparison of 
expansion of G{z,z') obtained from eqs. (^), (^), and similar equation for 1^2, in the region to 
the left of the junction with the expansion to the right establishes four relations between the 
amplitudes E and F in these two regions: 

t'\F+{+0) - E+{+Oy) = (F+(-0) -rti?_(-0))t, 
(F_(+0) - r'^E+{+0))t' = t^(F_(-0) - E_{-0)r) , (5) 
t'fE^{-0)t = E^{+0) + r'^E+{+Oy - F_(+0)r' - r'^F+{+0) , 
t'^E+{+0)t' = E+{-0) + r^E_(-0)r - F+(-0)r - r1'F_(-0) . 

"Rotating" the amplitudes E and F by the unitary matrices uj we can bring all coefficients in 
eqs. (^) into the diagonal form in the transverse mode space. The equilibrium Green's functions 
Gi^2 deep inside first and second superconductors, which play the role of the source terms in 
equations defining the amplitudes, are proportional to unit matrix in this space. This means 
that the boundary conditions @ become diagonal for the appropriately rotated amplitudes. In 
terms of the combinations of these amplitudes 

Qf = {u4F+i-0)u\±uiF^i-0)ul))/2, Qf = (n3F+(+0)4 ± ^^2i^-(+0)4))/2 , 
Pf = (u4^+(-0)4 - uiE^{-0)ul)/2 , P+ = (n4^+(-0)4 + uiE^{-0)ul)/2 - Gi , 
P2" = (u3^+(+0)4-n2^-(+0)4)/2, P+ = (u3^+(+0)4 + ^^2^-(+0)4)/2-G2, 

the boundary conditions @ can be written as follows: 

Pi = P2 , Qi = Q2 ; 
Q+-Q+ = + P+ + Gi + G2) , (6) 

G1-G2 + Pi+ - P^ = v^iQt + Qt) . 

Requirement that all deviations of the Green's functions from the equilibrium homogeneous 
values Gi^2, that are represented by the non- vanishing amplitudes P and Q, decay away from 
the junction region (at z — > ±00) imposes another conditions on P and Q. This condition has 
the same form as in the single- mode case 13: 

GjPj^ = -Pj^G, = {-iy^'P^, GjQf = QfG, = {-iy+'Qj, j = l,2. (7) 



Combined, eqs. (^) and (|^) give a system of linear matrix equations for the amplitudes P and 
Q. Solving this system and using the normalization condition satisfied by the Green's functions 
^1,2 ; ^1,2 = g^t finally: 

E+ = u\D{Gi + l)G-G+/Kui + Gi, F+ = u\^{Gi + I)/ Km , (8) 
E_ = u\d{Gi - l)G-G+/Kui + Gi , F_ = u\^/R{Gl - l)/Kui , 

where G± = (Gi it G2)/2 and K = 1 — DG^. The amplitudes E and F in eq. (|8[) are given in 
the first electrode {z — > —0). Similar expressions can be obtained for the second electrode. 

Equations (^) solve the problem of electron transport in different types of junctions ranging 
from quantum point contacts with few propagating electron modes to short multi-mode SNS 
junctions with diffusive electron transport. They are valid for arbitrary ratio of the two values 
of the superconducting energy gap in the junctions electrodes, and can be applied, e.g., to 
symmetric SNS junctions with the same gap and NS junctions with vanishing gap in one of 
the electrodes. The amplitudes (jsl) define the Green's functions (|2|) in the junction region and 
allow one to calculate different transport characteristics of the junction. For example, since the 
single-mode version of eqs. (|8|) was the starting point of the theory! of ac Josephson effect in 
high-transparency junctions, eqs. (|8|) prove that this theory can be extended automatically to 
the multi-mode junctions. Another application of eqs. (^) is presented in the next Section, where 
they are used to calculate equilibrium fiuctuations of the supercurrent. 



3 Equilibrium supercurrent noise 

General expression for the current noise in terms of the quasiclassical Green's functions is derived 
in 1121. In equilibriurn, this equation reduces to the following result for the spectral density of 
current fiuctuations El: 

Siiu)) = 32^2;^ J def{e){l- f{e±uj))Tr[p±{e)azP±{e±uj)az-u±ie)aziy:f{e±u;)az], (9) 

where az is the Pauli matrix, the trace is taken over the transverse mode space and over electron- 
hole space, p(e) = E^{e) — E^{e), and z^(e) = F^{e) — F^{e). The functions with superscripts 
R, A denote retarded and advanced components of the Green's functions. 

In what follows, the noise (jH) is calculated for a symmetric junction between two identical 
superconductors which are close to "ideal" BCS superconductors: 



Here (p is the Josephson phase difference across the junction. Imaginary parts 71^2 of e, A arise 
from inelastic scattering, for instance the electron-phonon scattering, and are assumed to be 
small, 71^2 ^ e , A. The functions G2 ' in the second electrode are given by the same expression 
with Lf — > —Lf. For superconductors described by the Green's functions (lTo|) , the amplitudes @ 
contain resonant denominators 1 — DG^, the contribution from which dominates the noise at 
low frequencies to <C A/h. Indeed, in this case G^ = A^ sin{(p/2) /(A^ — e^), where A = A it 172 
and e = e lb i^i, and 1 — DG^ nearly vanishes for e ~ = ibA[l — Dj^ sin^((/?/2)]^/^. Then, 

p(^),Ke)ocE(^^J;f^^2/4> (11) 



where 7^ = 2(7i(efc) =F72(eA;)efc/A). (Since 7i(-e) = 71(e) and 72(-e) = -72(e), the width 7 of 
the resonance is the same for both resonances in eq. (|Tl|).) The two resonances correspond to the 



two discrete Andreev states per mode with energies zte^ in the subgap region. Their broadening 
7 is caused by the inelastic transitions between these states and the continuum of states in tl^e 
bulk electrodes. Expression for 7 in the case of electron-phonon transitions is given, e.g., int 



Due to resonance (11), we can take all non-resonant terms out of the integral over e in eq. 



(|^. Evaluating the non-resonant terms at e = ite^. and intergrating the resonant denominators 
we obtain the spectral density of current fluctuations at low frequencies w ~ 7: 

1 ^ / T \ 2 



^^^'^) = 2^gUosh(.,/2T)j Z;^' (^2) 

where 1^ = (eA^/2?i)L'fc sin (/j/efc is the contribution to the supercurrent from one of the Andreev 
states in the kih. mode. 



Equation (12) determines the equilibrium noise of the supercurrent in short Josephson junc- 
tions. It has a simple interpretation in terms of the fluctuations of the occupation factors of 
Andreev states carrying the supercurrent. All states are occupied independently one from an- 
other. If a state with energy itefc is occupied, it gives contribution ±1^, to the supercurrent; if it 
is empty, the contribution to the supercurrent is zero. These two situations are realized with the 
probabilities f{±ek) and 1 — /(zbe^), respectively, and the transitions between them occur with 
the characteristic rate 7. Calculation of the spectral density of this simple classical stochastic 
process reproduces eq. (^). This equation was conjectured inEJ on the basis of the similar inter- 
pretation of the noise in the single- mode ballistic junctions. The calculation presented in this 
work provides its rigorous proof. 

At temperatures close to the superconducting critical temperature Tc of the junction elec- 
trodes, the energy gap is small, T ^ A, and the transition rates 7^ and the occupation prob- 
abilities of Andreev states are independent of the state energy, and eq. ( |T2[ ) can be simplified 
further. In this high-temperature regime, one also has to take into account the non-resonant 
contribution to the current noise from the continuum part of the junction spectrum which can 
become comparable to the resonant contribution (^). As can be seen from eq. (P), this non- 
resonant contribution coincides at T S> A with the regular equilibrium thermal noise of a normal 
junction. Combining the two contributions we obtain the noise in the high-temperature limit: 



5K-) = ^(^L4^j^7^ + ^^> (13) 

where Sn = GT/vr and G = (e^/vr^) J2k=i is the junction conductance. 

Equation (|T3| ) can be evaluated if the junction has known distribution of transparencies. An 
important example is provided by the regular multi-mode disordered SNS junction with diffusive 
electron transport. This type of junctions is characterized by the quasicontinuous Dorokhov's 
distribution of transparencies, J2k = (vr?iG/2e^) Jq dD/ D{1 — D)^/^. Calculating the integral 
over D we find spectral density of current noise at large temperatures: 

Siiu) = ^ cos2 ^ - 1) + Sn, ^ E [-vr, vr] . (14) 

2n 2 \smip / tj^ -|- 7^ 

The phase dependence of the supercurrent-noise part of eq. (|T^ ) is plotted in Fig. 1. Equa- 
tion (|1|) shows that if T <C A^/hj, this noise gives the dominant contribution to the cur- 
rent noise at frequencies tu The noise can also be larger than the average supercurrent 
I{f) = (7rGA^/4eT) sin If A/T < (e"^ /hG)^^"^ , typical r.m.s. of the supercurrent noise (in the 
frequency range 7) is larger than I{(p)- This means that the supercurrent noise should qualita- 
tively change the behavior of the SNS junctions at temperatures close to Tc and should also be 
observable at low temperatures. 
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Figure 1: Dependence of the zero-frequency spectral density S{0) — {Si{0) — Sn)/{GA^ /2h'y) of the supercurrent 
noise in short disordered SNS junctions on the Josephson phase difference (p at large temperatures T 2> A. All 

quantities are defined as in the eq. (|lj). 

4 Conclusion 

In conclusion, the equilibrium supercurrent noise has been studied in short Josephson junctions 
with arbitrary electron transmission properties. Characteristic qualitative features of the noise 
are its dependence on the Josephson phase difference across the junction, and the cutoff fre- 
quency determined by the relaxation rate, which can be much smaller than the temperature and 
superconducting energy gap. The r.m.s. of the noise is larger that the average supercurrent at 
temperatures close to the critical temperature of the junction electrodes. A formal result of this 
work is the general expression for the quasiclassical Green's functions of the multi-mode junc- 
tions that can be used to study other transport characteristics, for instance, the non-equilibrium 
noise in the MAR regime. 
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